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We investigate the effects of matter-enhanced (MSW) transformation of neutrinos in the Early 
Universe on the primordial helium yield {Yp). We find that Yp is affected much more by the MSW- 
induced alterations in the neutrino energy spectra than by the associated change in expansion rate. 
Specifically, the absence due to transformation of low energy electron neutrinos can significantly 
affect neutron-proton weak interconversion rates through the lifting of Fermi-blocking of neutron 
decay at low energies and through halting low-energy neutrino capture on neutrons. We find that 
the change of Yp within a causal horizon is —0.005 < ^V^ < 0.013 for — m^^ < 10* eV^ in the 

case of Vfj,^r-fs-Ve mixing, with the lower limit at m^^ ^ ^^^^s ~ 100 eV^, and —0.002 < SYp < 0.020 
for m^_, — m^^ < 1 eV^ in the case of Ve-Va mixing. 



I. INTRODUCTION 
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Big Bang Nucleosynthesis (BBN) remains one of the most successful probes of early times in the hot big bang 
' cosmology. Standard BBN (SBBN) assumes the Standard Model proposition of only three massless light neutrinos, 
^■f^ ' leaving the primordial element abundances characterized by only one parameter, the baryon-to-photon ratio 77. Now, 
the evidence for neutrino mass from the Super-Kamiokande atmospheric neutrino observations is overwhelming |Q] . In 
addition, the solar neutrino deficit and the LSND experiment |2| are hinting at the presence of a fourth light neutrino 
' that would necessarily be "sterile" due to the Z decay width [ 3 1 . The role of neutrino masses and their mixing bring 
■ another aspect to the physical evolution of the Early Universe, when neutrinos played a large role. The effects of 
neutrinos on big bang nucleosynthesis can be in part parameterized by the effective number of neutrinos (i.e. the 
, expansion rate), and in the weak physics which determine the neutron-to-proton ratio. In many works is regarded 
as the sole determinant of neutrino effects in BBN. However, we find that the effect of neutrino mixing on BBN is not 
well described by this one parameter. Characterizing the effects of neutrino mixing on BBN by Ni, can become overly 
" rS i' complicated, since the parameter has been used to describe several distinct physical effects, including lepton number 
Q , asymmetry and change in energy density. Another feature we find is that matter-enhanced mixing of neutrinos in the 
^ Early Universe can alter u,, spectra to be non-thermal, which uniquely affects the nucleon weak rates by lifting Fermi 
^ ' blocking of neutron decay by neutrinos and suppression of neutrino capture on neutrons. 

The observation of the abundances of primordial elements has greatly increased in precision within the past few 



years In SBBN, each observed primordial abundance corresponds to a value of the baryon-to-photon ratio 

' r] ss 2.79 X 10^r2^^/i~^, where fib is the fractional contribution of baryon rest mass to the closure density and h is 
■ the Hubble parameter in units of 100km s^^ Mpc^^. The theory is apparently successful since the inferred primordial 
abundances correspond within errors to a single value of 77. SBBN took a new turn in 1996 when it was found 
that the deuterium abundance at high redshifts may be significantly lower than that inferred through the cosmic 
helium abundance |^ . The advent of high-resolution spectroscopy of quasars allowed the measurement of deuterium 
abundance in high-redshift clouds between us and distant quasars. The "early days" of these high-redshift spectral 
deuterium measurements were filled with controversy due to a disparity in the inferred deuterium to hydrogen (D/H) 
ratio between two observational groups working with different quasars. More recently, the number of quasars with 
low D/H has grown. The deuterium abundance has converged to ~10% precision: (D/H)p ~ (3.4 ±0.25) x 10~^. The 
•q value inferred from this observation is (5.1 ± 0.5) x 10"^". The prediction for the corresponding primordial 

■^He mass fraction is Yp{D) ~ 0.246 ± 0.0014. 

Based on stellar absorbtion features in metal-poor HII regions, the observed helium abundance has been argued by 
Olive, Steigman, and Skillman (OSS) § to be Yp{OSS) ~ 0.234±0.002(stat) ±0.005(sys). The corresponding baryon- 
to-photon ratio is rjoss — (2-1 ± 0.6) x 10"^". The discrepancy between the central values of Yp{D) predicted by the 
Buries and Tytler deuterium observations and Yp{OSS) has led to much speculation on non-standard BBN scenarios 
which could reconcile the two measurements by decreasing the Yp{D) predicted by BBN. However, the systematic 
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and statistical uncertainty in the Yp{D) and Yp{OSS) puts them weh within their 2(t errors. Furthermore, Izotov and 
Thuan |lO) exclude a few potentially tainted HII regions used in OSS from their analysis and arrive at a significantly 
higher Yp ~ 0.244 ± 0.002(stat). The systematic uncertainty claimed in observed Yp may also be underestimated 

It has been suggested that this dubious Yp and D/H discrepancy with standard BBN could be hinting at the 
presence of sterile neutrino mixing in the Early Universe However, we find that the overall effect of active-sterile 
mixing in the Early Universe would only increase Yp and thus increase the possible disparity. 

Active-sterile mixing can produce an electron lepton number asymmetry L{ve) before and during BBN p^ , p^ . If 
the sign of L{ve) is the same throughout the Universe, the helium abundance would change to be more favorable 
{SYp < 0) or less favorable {SYp > 0) with observation. Taking into account the effects of energy density, lepton 
number generation and alteration of neutrino spectra, we find that SYp can be large in the positive direction, but 
limited in the negative direction fl^ . In particular, we find that matter-enhanced transformations can leave non- 
thermal spectra which significantly affect the weak rates through the alleviation of Fermi-blocking and suppression 
of neutrino capture on neutrons Furthermore, causality limits the size of the region with a certain sign of L[ve) 
to be within the size of the horizon during these early times (~ 10^" cm at weak freeze-out) The sign of the 

asymmetry will vary between regions, and the net effect of transformation on Yp will then be the average of both the 
L{ve) > and L{iye) < BBN yields, which leaves SYp > 0. 



Two possible schemes exist for producing a nonzero L{Uf.) involving matter-enhanced (MSW) transformation of 
neutrinos with sterile neutrinos. These two schemes differ greatly in the energy distribution of the asymmetry in the 
spectra of the v^jvf, the asymmetry. The resonant transformation of one neutrino flavor to another depends on the 
energy of the neutrino. In the case of BBN, the active neutrinos all start off as thermal, and sterile neutrinos are 
not present. As the Universe expands, decreases in density and cools, the position of the MSW resonance moves in 
energy space from lower to higher neutrino energies. The energy of the resonance is |T4] 



where (£^/T)ros is the neutrino energy normalized by the ambient temperature T. The asymmetry L(va) is generated 
as the resonance energy moves up the neutrino spectrum. As the Universe cools, the resonance energy increases. 
This produces a distortion of the neutrino or anti-neutrino spectrum. There are two plausible cases for producing an 
asymmetry L(y^) that would affect the production of primordial Helium through the weak rates: 

• Direct Two Neutrino Mixing: An asymmetry in L(y^ is created directly through ^ Vf^ ^ Vg or v,, Vg 
resonance. In this case, the resonance starts at low temperatures for low neutrino energies, when neutrino- 
scattering processes are very slow at re-thermalization. This case leaves the i^e or spectrum distorted, with 
the spectral distortion from a thermal Fermi-Dirac form residing in the low energy portion of the spectrum. 
The position of the distortion cut-off moves through the spectrum before and during BBN. In calculating the 
effects of this scenario on BBN, we included the evolving non-thermal nature of the neutrino spectrum. 

• Indirect Three Neutrino Mixing: An asymmetry in Liv^) or L(yp) created by a z^r ^ Vgiv^;^ — > z^s) 
resonance is later transferred into the L(y^ through a — > v^(y^^ — > resonance (L(y^ > 0). This will 
happen for the anti-neutrino flavors for the Liy,,) < case. In this scenario, the resonant transitions occur at 
higher temperatures when neutrino-scattering is more efficient. Both the z^r('^/i) and z/e spectra re-thermalize, 
but since the temperature is below weak decoupling, the asymmetry in neutrino lepton numbers remain. The 
effect of this scenario on BBN can be described by suppression or enhancement over the entire Ve/i^e spectra. 



II. MATTER-ENHANCED MIXING OF NEUTRINOS AND NEUTRINO SPECTRA 

IN THE EARLY UNIVERSE 
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III. NEUTRON-PROTON INTERCONVERSION RATES AND NEUTRINO SPECTRA 



BBN can be approximated as the freeze-out of nuclides from nuclear statistical equilibrium in an "expanding box." 
Nuclear reactions freeze-out, or fail to convert one nuclide to another when their rate falls below the expansion rate 
of the box (the Hubble expansion). As different reactions rates freeze-out at different temperatures, the abundances 
of various nuclides get altered and then become fixed. The nuclide produced in greatest abundance (other than 
hydrogen) is *He. The reactions that affect the ^He abundance (Yp) the most are the weak nucleon interconversion 
reactions 

n + v^, ^ p + e~ n + p + n p + e~ + (2) 

These reactions all depend heavily on the density of the electron- type neutrinos and their energy distributions . 
The two neutrino transformation scenarios for altering the i>e spectrum described above affect Yp through these rates 
and through the effect of increased energy density on the expansion rate. 



A. Direct Two-Neutrino Mixing: Vs «-> Vs 



In this situation, cither the populations of the fe or i/g energy spectra are cutoff at lower energies. The reaction 
rates that alter the overall n ^ p rates the most because of the spectral cutoff are 



= aJ v,ElEldp,[l + e-^^^A^^]-i[i + ^-EJkT.^-i (3) 



\n.^,e =Aj v^EyMe"""^"^-^ + ir'[l + e-^=/^-^-]-^ (4) 

(the notation in these expressions follows that in Ref. jlj]). The rate in Eq. || is modified with the Pe spectra, and 
the rate in Eq. ^ is modified with the spectra. 

For a spectral cutoff for z7e, corresponding to Ds ^ f'e transformations, the Fermi-blocking term [1 + ^-^.y/kT^^-i jjj 
neutron decay (Eq. ^ becomes unity at low P^, energies. The rate integrand is then greatly enhanced when low energy 
neutrinos would otherwise block the process (see Figure |l|) . In the case of a spectral cutoff for i/f. , corresponding to 
i^s ^ V(, transformations, the spectral term p1[e^"^^'^'^ + 1]^^ and thus the rate integrand go to zero for low energy 
neutrinos (see Fig|3). 
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FIG. 1. Plotted are the contours of the amplitude of the rate integrand for n p -f e~ + i^e, with increasing amplitude 
towards the center, (a) shows the lack of fermi blocking by low energy electron anti-neutrinos at higher temperatures, (b) 
shows the standard BBN integrand. 
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FIG. 2. Plotted are the contours of the amplitude of the rate integrand for n + p + , with increasing amplitude 

towards the upper right: (a) shows the blocking of the integrand for lower energies due to the suppression of low energy electron 
neutrinos; (b) shows the standard BBN integrand. 

The resultant change in Yp for this "two neutrino" case is shown in Figure ^. For De Dg transformations, 
L(ve) > 0, the neutron decay rate in Eq. |^ is enhanced, fewer neutrons are present, and Yp is decreased. For Vs 
transformations, i(t'e) < 0, the n z/g — > pe rate in Eq. ^ is suppressed, more neutrons are present, and Yp is increased. 
The effect of the rate in Eq. |4| is much greater than that of the rate in Eq. || on the neutron-to-proton ratio and 
Yp. This causes a larger magnitude change in the neutron-to-proton ratio for the L{i'e) < case, and thus a larger 
magnitude change in Yp in the positive direction than in the negative. 

The sign of the asymmetry Liy^) is randomly determined [p^ . Because of causality, different horizons will have an 
asymmetry of random sign fifty-percent positive and fifty-percent negative. The effect of the asymmetry on Yp 
will then be the average of these two cases, which is also plotted in Figure ||. 
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FIG. 3. Plotted are 5Yp vs. 5m^ for the direct two-neutrino transformation scenario. The top line is for L(ve) < 0; the 
bottom is for Liue) > 0; the middle line is the average of both cases enforced by causality. 

B. Indirect Three-Neutrino Mixing 

In the case that an asymmetry is produced in Ve/ve indirectly through an initial transformation between Vt 
or — > Ug, the neutrinos can re-thermalize through self-scattering at higher temperatures. The number densities 
remain asymmetric, but the asymmetry is spread uniformly over the entire energy spectra. The effect on the weak 
rates is not as drastic as the above case, but is still present. In this case, the increase in energy density due to the 
thermalization of sterile neutrinos also affects Yp through the increase in the expansion rate. Thus, for L{ve) > 0, 




the suppression of Yp is counteracted by the increased expansion rate. The increase in energy density always causes 
a positive change in 1^, while the change caused by the weak nucleon rates depends on the sign of L{i>e), thus the 
overall resulting change in Yp is greater in the positive direction than in the negative (see Figure^). For the L(i^e) > 
case, the negative change SYp has a minimum at Sm^ « 100 eV^ since at this point the positive effect of energy density 
begins to dominate the negative effect of the modification of the weak rates. 
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FIG. 4. Plotted are SYp vs. 5m^ for the indirect three-neutrino transformation scenario. The top line is for L{ue) < 0; the 
bottom is for L{ve) > 0; the middle line is the average of both cases enforced by causality. 



IV. CONCLUSIONS 

The characterization of the effects of neutrino transformations on BBN through a single parameter Njy hides varied 
physical effects of neutrino mass and mixing that are non-trivially related. These effects include an increase in 
energy density and modification in the nucleon weak rates. The modification of these rates is through either an 
asymmetry that is distributed over the entire spectrum, or an asymmetry confined to low energy v^jv^. We find 
that the distribution of the L(ye) in the spectra is important to the weak rates, the neutron-to-proton ratio, and the 
production of *He. In addition, the change in the predicted Yp is always greater in the positive direction than in the 
negative, thus when averaged over causal horizons, SYp is always positive. The changes SYp are comparable to the 
uncertainty in current Yp measurements, thus resonant neutrino mixing can measurably affect the BBN predictions. It 
is important to note that there may remain non-trivial portions of parameter space in non-standard BBN — including, 
but not limited to, neutrino mixing, spatial variations in 77, and massive decaying particles — that fit within the 
observed uncertainties of the primordial abundances. 

K. A., X. S. and G. M. F. are partially supported by NSF grant PHY98-00980 at UCSD. 



[1] Y. Fukuda, et al, Phys. Rev. Lett. 81, 1562 (1998). 
[2] D. O. Caldwell, Int. J. Mod. Phys. A13, 4409 (1998). 

[3] Particle Data Group, Review of Particle Physics, European Physical Journal C3, 1 (1998). 



[4] D. Schramm and M. Turner, Rev. Mod. Phys. 70, 303 (1998); S. Buries et a/., |astro-ph/9901157 
[5] S. Buries and D. Tytler, Astrophys. J. 499, 699 (1998). 
[6] S. Buries and D. Tytler, Astrophys. J. 507, 732 (1998). 
[7] K. Olive, G. Steigman and E. Skillman, Astrophys. J. 483, 788 (1997). 
[8] P. Bonifacio and P. Molaro, Mon. Not. R. Astron. Soc. 285, 847 (1997). 
[9] D. Tytler, X. Fan and S. Buries, Nature 381, 207 (1996). 
[10] Y. Izotov and T. Thuan, Astrophs. J. 500, 188 (1998); ibtd. 497, 177 (1998). 



[11] G. J. Mathews, R. N. Boyd and G. M. Fuller, Astrophys. J. 403, 65 (1993). 
[12] R. Foot and R. R. Volkas, Phys. Rev. D 56, 6653 (1997). 
[13] X. Shi, Phys. Rev. D 54, 2753 (1996). 

[14] X. Shi, G. M. Fuller and K. Abazajian, Phys. Rev. D 59, 063006 (1999). 
[15] K. Abazajian, X. Shi and G. M. Fuller, in preparation. 

[16] X. Shi and G. M. Fuller, submitted to Phys. Rev. D. 

[17] R. A. Alpher, J. W. FoUin, Jr., and R. C. Herman, Phys. Rev. 92, 1347 (1953). 
[18] S. Weinberg, Gravitation and Cosmology, J. Wiley, New York, 1972. 



6 



